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We suggest a mechanism which promotes the existence of a phase soliton - topological defect 

formed in the relative phase of superconducting gaps of a two-band superconductor with s^ type 

of pairing. This mechanism exploits the proximity effect with a conventional s-wave superconductor 

which favors the alignment of the phases of the two-band superconductor which, in the case of s.\ 

pairing, are 7r-shifted in the absence of proximity. In the case of a strong proximity such effect can 
be used to reduce soliton's energy below the energy of a soliton-free state thus making the soliton 
thermodynamically stable. Based on this observation we consider an experimental setup, applicable 

both for stable and metastable solitons, which can be used to distinguish between s^ and s++ 

types of pairing in the iron-based multiband superconductors. 



Introduction. In the last decade multiband supercon- 
ductivity became a central topic in condensed matter 
physics. The discovery of MgB2 (Ref. [Tj) and, more re- 
cently, of an entire family of iron-based high-temperature 
superconductors (Ref. EJ) showed that it is not an exotic 
possibility, but a problem with enormous theoretical and 
practical relevance. Although the origin of superconduc- 
tivity in MgB2 and iron pnictides appears to be quite dif- 
ferent (phonon-driven in the first, and purely electronic 
in the case of the latter), in all of these materials the pres- 
ence of multiple gaps leads to important and far-reaching 
consequences. In fact, it has been argued that the multi- 
band character of the Fermi surface is essential for the 
superconductivity in pnictides. 

It has been long known that two-band superconductors 
(SCs) have excitations associated with the fluctuations 
of the relative phase of the two gaps known as Leggett 
mode[3j. Only recently, however, has it been recognized 
that there are also topological excitations associated with 
the multiple gaps[4j-(9]. In their simplest version these 
excitations can be thought of as soliton-like domain walls 
between regions in which the relative phase differs by 2n. 

In this paper we argue that such solitons can be used to 
distinguish between the conventional s ++ and the more 

exotic s-| pairing in which phases on the bands are 7r- 

shifted. Both states have been suggested as a possibility 
for iron pnictides. Despite being physically very different, 
distinguishing them experimentally is not a trivial task 
since they belong to the same symmetry class (for some 
interesting suggestions see e.g. Refs. fT0l-H2]l. Here we 
elaborate a practical setup which can be used to detect 
the difference between the two order parameters. It uti- 
lizes the proximity effect which affects the existence of a 
soliton in a qualitatively different way depending on the 
relative sign between the gaps. 

We consider a ring made of a two-band SC, segment 
of which is covered with a stronger, higher T c , s-wave 



superconductor ("proximity patch"). By the proximity 
effect the phases of the two bands will tend to align with 
the phase of the s-wave side (as discussed later) , which 
reduces the energy of the soliton if the two-band SC is in 

s_| state. For sufficiently strong proximity, the energy 

of the soliton can be brought below the energy of a vor- 
tex, thus making the former thermodynamically stable. 
We suggest to look for both equilibrium and metastable 
solitons by using magnetization measurements. 

Description of the model. Phenomenological descrip- 
tion of a two-band superconductivity requires a two- 
component order parameter [TH - HB] . Let m, be the effec- 
tive mass of charge carriers in the i-th band and A denote 
the vector potential. Following early works we consider 
a two-component Ginzburg-Landau (GL) model [17] 
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In this model each of the bands is assigned a com- 
plex gap ipi = \vbi\e lBi and is described by a standard 
GL functional with the potential energy term l^(lV'il) 
(first line in Eq. (flj)). Coupling between the bands 
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FIG. f : Proposed experimental setup which utilizes a prox- 
imity patch for the observation of a phase soliton. 



is provided by a Josephson energy Ej which depends 
on their relative phase (second line in Eq. (fil). In 
terms of the band phases this coupling can be written 
as Ej = J —2K\tpi\\ip2\ cos(#i — 9 2 ). Sign of the coupling 
constant k determines the relative phase of the bands in 
the ground state: k > favors 0\ = 02 i.e. the phases 
are aligned while k < favors anti-aligned configuration 
Oi = 02 + 7T. These states are often referred to as s ++ 
and s_| states, respectively. 

Generally speaking the possibility of the relative phase 
dynamics depends on the strength of n. If |k| is very large 
then the phases are locked to each other and, as a result, 
act synchronously. However for weaker \k\ behavior of 
the phases can be different. In particular one can envision 
a state in which phases of the order parameters on the 
two bands wind differently e.g. host different number of 
vortices. A conceptually similar state - half-quantum 
vortex - has been recently observed in spin-triplet SC 
Si'2Ru04 where the role of the bands is played by the 
two spin components of charge carriers [15]. 

To further elaborate on this idea we consider a ring of 
radius R, threaded by magnetic flux $ (see Fig. HI . The 
height w and the thickness d of the ring should be compa- 
rable or smaller than the characteristic lengthscale of the 
spatial variations of the order parameter [19]. The single- 
valuedness condition requires both phases 0i to wind only 
modulo 2-7T when taken around the ring. For the relative 
phase a = 9± — 82 this implies that 



a(x + 2ttR) = a(x) + 2im ai 



(2) 



where x is the coordinate along the ring. For a regu- 
lar vortex state n a = and the two phases have the 
same windings 2D]. Importantly, there also exist states 
for which windings of #1,2 around the ring are different 
corresponding to n a — ±1,±2... and non-zero wind- 
ing in a. Competition between the kinetic energy which 
tends to spread such winding and the Josephson cou- 
pling which, independent on its sign, tends to localize it, 
leads to the appearance of a spatially confined profile - 
kink - of the gradient of a. Such solution is called phase 
soliton [H [5T]. Our aim is to show that the stability of 
the phase soliton can be affected by the proximity effect. 
Proximity effects between conventional s-wave and 

s_i two-band SCs have been studied recently [2"2Tl24| . 

While in the case of proximity between conventional 
(even multigap) SCs the effects are rather trivial, the 

proximity with s_| order parameter provides for much 

richer physics. There is a number of anomalous features 
which can appear in such structures. In particular, it has 
been found that under some circumstances the phases on 
the two bands tend to align with the phase of the s-wave 
SC. As a results the relative phase between the bands 
deviates from the bulk 7r value (note that the resultant 
state is complex and thus breaks time-reversal symme- 
try). If the proximity effect is sufficiently strong it is 
even possible to reduce the phase difference to zero and 



induce an s++ state in the (originally) s_| system. 

It is intuitively clear that the proximity effect described 
above can affect soliton's stability. Indeed, in the case of 

s_| symmetry the deviation of the relative phase from 7r 

(configuration < >) - essential requirement for the exis- 
tence of the soliton - is promoted in the region of the 
contact, thus lowering soliton's energy. It is also clear 
that if the contact is made between s ++ and s-wave ma- 
terials this will, at best, preserve soliton's energy since 
any deviation from 'f j" configuration is now further pe- 
nalized in the region of the contact. 

The qualitative discussion of the soliton's energetics 
given above is quite general and relies only on the exis- 
tence of the soliton and the proximity effect. To illustrate 
our arguments quantitatively we turn to a simplified, an- 
alytically solvable model, based on the framework of the 
GL theory ([l]). In it the effect of the patch is mimicked 
by a spatially dependent interband Josephson coupling 
k{x). The profile of this effective n(x) along the ring 
should be chosen to promote the alignment of the phases 
in the region of the patch. More precisely, we require 



^in ^ ^out 



(3) 



where «; n = k(x) for x in the patch and K out = k(x) 
otherwise. Notice that while for an s ++ ring k(x) has 

the same (positive) sign for all x, for an s_| ring where 

Kout < 0, k(x) can, in the case of the strong proximity, 
reverse its sign in the patch area so that KoutKm < 0. 

Before starting on the calculations we should discuss a 
potential caveat related to the use of the GL free en- 
ergy expansion. It has recently been argued [23] that 
the accuracy of expansion (HI which is controlled by the 
reduced temperature (T c — T)/T c imposes a constraint 
tpi = real number xip 2 (see, however, Ref. |26, ). All other 
solutions, including the soliton one, although not forbid- 
den in principle, should be obtained by supplementing 
expansion |l]) with higher-order terms [2"T] (cf. Ref. 128)). 
We should point out that such argument is inapplicable 
in the presence of the proximity patch - the main in- 
gredient of our proposal. Indeed, it is quite clear that 
independent of a particular form of description the rela- 
tive phase far from the patch is that of a bulk material 
while deep in the region of the patch a is controlled by the 
proximity effect as shown by a microscopic analysis [2^]. 
The existence of a (meta)stable soliton in our model re- 
lies only on an interpolation between these limits and is 
thus highly plausible. 

To continue our analysis we assume that the magni- 
tudes of the order parameter \ipi\ is field- and coordinate- 
independent which allows one to drop the potential en- 
ergy V^(|-0i|). Let us introduce the following notation: 



-2 _ \-2 



X7' = 87re^i|7Kc^), A~ 2 = A^ + A 



(4) 



We focus on the limit of weak screening defined by 
Rd/X 2 <C 1. In this limit the difference between ap- 




FIG. 2: (a) Spatial profile of the relative phase of the soliton 
and non-soliton solutions for s+_ pairing in the presence of 
the proximity patch, (b) Dependence of the excessive soliton's 
energy on the length of the proximity patch for several K. 

plied and total fluxes, as well as the difference between 
the Gibbs potential and the free energy, can be ignored. 
Following Ref. [13] we define new phase variables 9 = 
\ 2 (9iX[ 2 + 02A2 2 ) and a = 0± — 6*2- It is now possible to 
integrate out the ^-dependence in Eq. pj which leads to 
the following expression for the free energy of the ring|29j: 



F = F* + F a (5) 

F* = eo [m(A 2 /A 2 ) + n 2 (A 2 /A 2 ) - $/$ ] 2 (6) 

fyd x a) 2 -K(x)cosa} (7) 



F a [a) = e ,„.„ / dx 



n\ 2 X 2 



where $0 = hc/2e, rii is an integer which describes the 
winding of 9i, €q = wd&o/(16n 2 R\ 2 ) is the electromag- 
netic energy scale and K{x) = k(x) mim 2 < I ) o/(27r 3 /i 4 A 4 ) 
is the renormalized Josephson coupling. Minimizing F a 
- the flux-independent contribution associated with the 
relative phase - we find that a(x) is determined by the 
sine-Gordon equation with a spatially dependent mass: 



d 2 a(x) — K{x) sin a(x) = 



(8) 



subject to the boundary conditions specified by Eq. (pi). 
At any given flux $ thermodynamically stable state is 
the one which realizes the global minimum of F. It can 
be shown from Eqs. ([5J)-(f7]) that in the absence of a prox- 
imity patch (i.e. when K(x) — Kq = const) the global 
minimum never corresponds to a soliton state. While -F$ 
itself allows for a stable soliton, its effect is countered by 
F a , which satisfies the following inequality: 



AF a ) > eo/4 



(9) 



where AF a = (F a [a] — F a [0])\ K=Ko is the difference 
between relative phase energies of a soliton and a con- 
ventional vortex solutions for K(x) — Kq. The lower 



bound for AFa is reached when Ai = A2 and K = 
i.e. for a uniform phase winding. In this limit electro- 
magnetic response of the soliton is equivalent to that of 
a half-quantum vortex [18 v . 

From now on we specialize on a physically reasonable 
setup in which lengths of both the patch and the soliton 
are small compared to the perimeter of the ring[30j. In 
this limit single-valuedness condition (pi) can be replaced 
with boundary conditions at infinity: a(oo) = a(— 00) + 
2ir for a soliton state and a(— 00) = a(oo) otherwise. We 
have excluded states with larger winding numbers as they 
have higher energy. In terms of the winding numbers of 
each band the above conditions imply |m — no\ < 1. 

We now demonstrate that constraint (pi), which pre- 
vents the thermodynamic stability of the soliton, can 
be lifted in the presence of the proximity patch. De- 
tailed calculation of K(x), describing the full effect of 
the patch, requires the use of the microscopic theory and 
is beyond the scope of this paper. Generally speaking 
function K (x) should be such that condition (pi) is sat- 
isfied. A straightforward analytical solution is available 
for a delta-patch model: 



K{x) = K (l - 7 |iY r 1/2 <5(x - 10)) 



(10) 



which describes a small, but very strongly coupled patch, 
parametrized by 7. For an s ++ ring we have Kq > and 

7 < 0, while for a non-trivial s_j pairing the opposite 

holds: K < and 7 > 0. 



The presence of the (5-function in Eqs. (10) and (18 



can be dealt with by imposing an additional boundary 
condition on the unperturbed (K(x) = Kq) solution ob- 
tained by integrating Eq. (pi) around the location of the 
(5-peak. We note that in this case a non-soliton solution, 
like a soliton one, may have a non-trivial structure in 
which d x a 7^ 0. Let us define AF a 1 ' to be the difference 
in F a evaluated between the lowest energy soliton and 
non-soliton solutions for the ring whose pairing symme- 
try type is i. A straightforward calculation shows that 



AF a ++^ = 8F a0 



APi-*~)=2F tt0 x{jr- ^" 2 M2, 



4/7, 



7 <0 

<7 

7 > 2 



(11) 



where F a0 = e i?| i^ | 1/2 A 4 /(7rAf A|). The above result is 
applicable when RIKq] 1 ^ 2 > 1. 



It follows from Eqs. (11), (12) (in agreement with the 



earlier discussion) that the effect of the proximity patch 
on the winding of the relative phase is qualitatively dif- 
ferent for ,s ++ and s^ pairing symmetries. While in 

the former case the a-winding of the soliton always costs 
non-zero energy AF a limited from below by Eq. (|9|, 

for s^ symmetry the energy difference AFa can be 

made arbitrary small, thus circumventing restriction (pi 
and rendering the soliton thermodynamically stable. The 
independence of AF++ on 7 is a consequence of the fact 




FIG. 3: Magnetic moment for stable, (a), and unstable, (b), 
soliton configurations in a ring geometry. Solid lines corre- 
spond to the equilibrium behavior. Dashed and dotted lines 
correspond to metastable soliton and vortex states respec- 
tively. 

the soliton's kink and the patch "repel" each other for 
s ++ pairing, and the lowest energy configuration corre- 
sponds to an infinite separation between them. 

The conclusions derived from the delta-patch model 
are confirmed by numerical calculations for a finite length 
proximity patch. In this model the proximity is induced 
in a finite segment of length L p such that K (x) = K' = 
const on the patch and K(x) = K = const otherwise. 
We focus on the limit when the length of the soliton's 
kink is much smaller than the perimeter of the ring. The 

results of the calculations for s^ pairing are presented 

on Fig. [2J Fig. [2^, shows the spatial profile of the relative 
phase for soliton and non-soliton solutions and Fig. [2)3 
demonstrates the reduction of the relative energy of the 
soliton as the length of the proximity patch is increased. 

Proposed experimental setup. We now come to the 
question of experimental observation of phase solitons. 
We consider the setup illustrated on Fig. [T] To insure 
strong effect of the proximity the thickness of the ring 
in the direction perpendicular to the interface should be 
comparable or smaller than the coherence length. 

There are several other factors which would favor 
strong proximity effect 131. • First, the s-wave SC has 
to have higher T c . Second, the interboundary coupling 
between all the gaps has to be sufficiently strong. This 
implies not only low boundary resistivity, but also re- 
quires the ratio of the normal state conductivities of the 

s and s-| materials to be much larger then one - thus 

ensuring that the effect of the s-wave SC on the 8-\ side 

is maximal [32 . Based on these requirement we suggest 
that a good candidate combination consists of a member 
of the 11 iron chalcogenide or 111 iron pnictidc families 
(which arc semimetals with T c 's of about 10K and typi- 
cally have gaps with close values; some of these materials 

has been suggested to be s_| , Refs. |33j [34]), coupled to 

MgB 2 as an s-wave material. MgB 2 has T c w 42K and 
is a good metal in its normal state. Its own multigap 
nature should not be a problem in this context (provided 



it is in s ++ state). 

In a ring geometry the phase soliton can be observed 
through measurements of ring's magnetic moment M. 
Fig. [3] shows the dependence of M on the applied flux <I> 
in the presence, (a), and in the absence, (b), of a ther- 
modynamically stable soliton, obtained from Eqs. (pi) - 
([7]). A stable soliton state is seen as two extra transitions 
between the adjacent vortex (fluxoid) states. These two 
transitions correspond to a mismatched phase winding 
in the two bands |m — 712] = 1 and are in general split 
because of the different band superfluid densities. 

Even if the coupling to the proximity patch is not suf- 
ficient to induce the thermodynamic stability of the soli- 
ton, it can still be observed as a metastable configuration. 
Metastability of vortices is well-known [35] and can be a 
notorious problem. We suggest to look for metastable 
solitons through thermal cycling, which in practice can 
be realized by heating a part of the ring with a laser 
pulse and then allowing the ring to cool down [3S1 [37J . 
To observe the soliton the system is repeatedly thermally 
cycled at a fixed value of the applied flux and a distribu- 
tion of magnetic moments which correspond to different 
metastable states is collected. In the absence of solitons 
such distribution plotted at different values of flux should 
form a set of equidistant lines (dotted lines on Fig. [31 see 
also Fig. 1 in Ref . l37j) . In the presence of a (metastable) 
soliton extra lines appear (dashed lines on Fig. [3]); the 
relative number of points on these lines gives the proba- 
bility of accessing a metastable soliton state |38j. 

The metastability of a non-equilibrium state can be 
greatly improved for a thin ring, where dynamics of the 
order parameter is one-dimensional. In realistic experi- 
mental settings the probability of accessing a metastable 
state with an excess energy E separated by the barrier 
V from the ground state is of the order of e~ xE l v where 
x is a parameter which depends logarithmically on the 
product of cooling time and characteristic attempt fre- 
quency. In the case of a thin ring interrupted by a weak 
link dominant contribution to E is of the electromagnetic 
origin so that E ~ $q/R 2 . At the same time the dom- 
inant contribution to V is set by the weak link's energy 
making V only weakly independent. Hence to observe 
metastable states the use of larger rings is desirable. 

Although the detection of the soliton for an s++ ring 
in the presence of the s-wave proximity patch is highly 
unlikely, it is nevertheless a possibility. To rule out such 
scenario the dependence of the stability region A<I> so i (for 
the thermodynamically stable soliton) or of the proba- 
bility of accessing a soliton through thermal cycling p so i 
(for a metastable soliton) should be checked for different 

lengths of the patch. In the s_| case the soliton's energy 

is reduced as the length is increased (Fig. [2b ) which leads 
to the growth of A$ so i and p so i, while for s++ it, at best, 
remains constant. The observation of such growth would 
unambiguously imply the s_| pairing. 

To conclude, we have considered the influence of the 



proximity effect on the stability of a phase soliton in a 
two-band superconductor. We have shown that the prox- 
imity between s^ and s-wave superconductors can re- 
duce the energy of the soliton and suggested a practical 
experimental setup which utilizes this observation to dis- 
tinguish between s^ and s ++ types of pairing. 
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